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We investigate non-Gaussianity in the modulated reheating scenario where nuc- 
lei ! tuations of the decay rate of the inflaton generate adiabatic perturbations, paying 
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particular attention to the non-linearity parameters /nLj^nl and (/nl as well as the 
scalar spectral index and tensor-to-scalar ratio which characterize the nature of the 
primordial power spectrum. We also take into account the pre-existing adiabatic 
perturbations produced from the inflaton fluctuations. It has been known that the 
non-linearity between the curvature perturbations and the fluctuations of the de- 
cay rate can yield non-Gaussianity at the level of /nl ~ C(l)> but we find that the 
non-linearity between the decay rate and the modulus field which determines the de- 
cay rate can generate much greater non-Gaussianity. We also discuss a consistency 
relation among non-linearity parameters which holds in the scenario and find that 
the modulated reheating yields a different one from that of the curvaton model. In 
particular, they both can yield a large positive /nl but with a different sign of $nl- 
This provides a possibility to discriminate these two competitive models by looking 
at the sign of #nl- Furthermore, we work on some concrete inflation models and 
investigate in what cases models predict the spectral index and the tensor-to-scalar 
ratio allowed by the current data while generating large non-Gaussianity, which may 
have many implications for model-buildings of the inflationary universe. 



1 Introduction 

Inflation is a promising candidate to generate primordial density fluctuations as well as 
to solve the horizon and the flatness problems. Although many observational supports 
for the inflationary scenario have been accumulated, it is still unclear what mechanism is 
really responsible for producing primordial density fluctuations. We usually assume they 
are generated by fluctuations of the inflaton, but those of other scalar fields can produce 
the primordial fluctuations too. Such scalar fields generically exist in the extensions of 
the standard model of particle physics, which motivates the curvaton scenario [1-3] and 
modulated reheating scenarios [4,5]. One of interesting features of these scenarios is that 
primordial non-Gaussianity can possibly be very large, in contrast to the case of single 
field inflation models where only the inflaton is responsible for density fluctuations so 
that almost perfect Gaussian fluctuations arise. Therefore, non-Gaussianity can be a very 
powerful tool to identify the real source of the primordial fluctuations. In fact, there 
has been reported that non-Gaussianity is detected in the cosmic microwave background 
almost at 3a level [6] although this is not confirmed by the latest WMAP 5- year results [7] . 
(See also Refs. [8-10].) Further observations of WMAP and future observations such as 
Planck can give us more information on non-Gaussianity and serve to discriminate different 
scenarios. 

In most works done thus far on this topic, only one source of the fluctuations has been 
considered. In general, however, since there exist a lot of scalar fields in a supergravity or 
superstring theory, different kinds of sources could contribute to the primordial fluctuations 
simultaneously. Therefore, it is interesting to consider models with mixed fluctuations from 
the inflaton and other sources. In such mixed scenarios, not only non-Gaussianity but also 
other features of the primordial fluctuations such as the amplitude, spectral index and 
tensor modes are affected in comparison to the case where the inflaton alone is the seed 
of fluctuations. Current cosmological data are already very precise regarding the scalar 
spectral index and the tensor-to-scalar ratio and they can severely constrain models of 
inflation. Some models of inflation are considered to have already been excluded [7]. For 
example, chaotic inflation with the higher order polynomial potential is excluded at more 
than 95 % confidence level. But, such constraints may be evaded by adding another source 
of fluctuations. In fact, it is shown that some models of inflation which are disfavored by 
the data can be liberated by adding the curvaton contribution [11-18]. In particular, in 
Ref. [17], it has been studied in what cases models of inflation can be relaxed by adding 
fluctuations from the curvaton in some detail assuming some concrete inflation models 
focusing on the spectral index and the tensor-to-scalar ratio. Furthermore, in Ref. [19], it 
was also discussed in what cases/models large non-Gaussianity can be generated satisfying 
the constraints on the scale dependence and tensor modes of primordial fluctuations in 
the mixed scenario. 

The modulated reheating scenario has been paid much attention recently as inter- 
esting other source of fluctuations and large non-Gaussianity but it is not so rigorously 
investigated as the curvaton mechanism in particular as regards mixed scenarios with the 



inflaton. Therefore, in this paper, we consider the mixed models where fluctuations from 
the inflaton and the modulated coupling can both contribute to the present cosmic den- 
sity fluctuations. Then, we study the effects of the contribution of the modulated coupling 
on inflationary parameters such as the spectral index and the tensor-to-scalar ratio and 
compare with the WMAP 5-year results. Furthermore, we also discuss in what cases large 
non-Gaussianity can be generated satisfying the observational constraints on the scale 
dependence and tensor modes of primordial fluctuations in such a mixed scenario. 

The structure of this paper is as follows. In the next section, we will give the expressions 
of the decay rate of the inflaton into radiation for various types of interactions between 
the inflaton and other fields and for various inflaton potentials, which will be used in the 
subsequent sections. In section [31 we will provide the e-folding number from the time when 
the current cosmological scales crossed the Hubble horizon during inflation to the time after 
the inflaton decays, which enables us to evaluate the curvature perturbations generated in 
this scenario. Then we give the expressions for the scalar spectral index, tensor-to-scalar 
ratio and three non-linearity parameters in section |U With the formalism summarized in 
section HI we work on some specific inflation models to compare the predictions of these 
inflationary parameters with observations, paying particular attention to in what cases 
non-Gaussianity can be very large. For inflation models which are considered to have 
already been excluded by the data, we also give discussions in what case the contribution 
from the modulated reheating can liberate the model. The final section is devoted to the 
conclusion and summary of this paper. 

2 Decay rate of the inflaton 

After inflation, the inflaton oscillates around the minimum of its potential and the period 
of the oscillations is much shorter than the expansion time characterized by the Hubble 
parameter. The energy density of the universe at that time is stored in the form of oscilla- 
tion energy of the inflaton. Since the universe should become radiation dominated before 
the time of big-bang nucleosynthesis, which is required to be consistent with observations, 
the inflaton must decay into radiation at some time after inflation. Regarding interactions 
between inflaton and radiation, we consider the following Lagrangian, 

£tat D " J2 VM^aA - Yl M a(^Wl ~ J2 ^(^Xl (1) 

a a a 

where <j) is the inflaton and Xa and ip a are scalar and spinor fields which constitute ra- 
diation (a represents the species of the particles). In the modulated reheating scenario, 
the coupling constants y a , M a and h a are functions of a scalar field which we denote as 
a. Although the number of such scalar fields is not necessarily one, we consider only one 
modulus to avoid inessential complexity. 

The oscillations of the inflaton act as a periodically changing external field on Xa and 
ijj a fields. This external field creates Xa or ip a particles out of the vacuum. Due to the 



total energy density conservation, the energy of created particles are compensated by the 
loss of the oscillation energy of the inflaton. Hence this process can be regarded as the 
particle production from the inflaton decay. 

Let us suppose that the inflaton potential around the minimum can be well approx- 
imated by a polynomial form as V(4>) oc <p 2n with n being a positive integer. Denoting 
the inflaton energy density as p^, the decay rate of the inflaton to the lowest order in the 
coupling constants is given by 

+ V ; ^ 8tt * ^ tormf *-" 8ir(mf) 3 
where rrf? is the effective mass of the inflaton defined by 



(ruff = V, 



(3) 



with (f> being the amplitude of the oscillations. Here and hereafter, the subscript <p of the 
potential V represents the derivative with respect to <p. A n , B n and C n are numerical 
coefficients of 0(1 ~ 100). The explicit values of these are given in the appendix 1X1 

3 The number of e-folding 

In the following, the e-folding number especially plays an important role in two aspects. 
First of all, when we calculate observables such as the amplitude, spectral index, non- 
Gaussianity and so on of the curvature perturbations, we have to know when the cosmo- 
logical scales exited the Hubble horizon during inflation, i.e. the e-folding number from the 
time of horizon crossing during inflation to the present. To this end, we need to know how 
the universe evolved from the inflationary universe into the standard radiation dominated 
universe, which obviously requires our knowledge of how the universe is reheated by the 
inflaton decay. Another aspect is that knowledge of the e-folding number enables us to 
calculate the curvature perturbations to any order in the perturbative expansion without 
invoking complicated perturbed equations but using the so-called SN formalism [20-23]. 
All the information regarding various observables of the curvature perturbations is con- 
tained in the e-folding number. In this section, we give an expression for the e-folding 
number from the time £* when the cosmological scale crossed the horizon during inflation 
to the time tf when the universe is completely reheated by the inflaton decay. After tf, 
we assume that the universe evolves according to the standard hot Big-Bang model. 

Let us denote the e-folding number from t* to tf as N(tf,t*,</>*,(T*). Here 0* and a* 
are values of and a at t = £*. For later convenience, we divide N(tf, £*, 0*, a*) into two 
parts as 

N(tf, U, 0*, <7*) = N(t end , £*, 4>*) + N(tf, tend, c*)j ( 4 ) 



where t en d is the time at the end of inflation. Hence, on the right-hand side, the first 
term represents the e-folding number during inflation and the second term remaining one 
after the inflation. In writing this equation, we have implicitly assumed that the mass 
and vacuum expectation value of the modulus are sufficiently small and the background 
dynamics during inflation is completely determined by the inflaton alone, which means 
that the first term on the right-hand side is a function of </>* only. Meanwhile, the second 
term depends only on er*. 

By using the slow roll approximation, we write the first term as 
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A(t end ,t*,0*) = / Hdt--—^ —d(j), (5) 

JU M P \ J 4,, v <f> 

where 0* and C nd represent the scalar field values at corresponding epochs. To give a more 
concrete expression for N(t cn d, £*, 0*), we need to specify the potential for the inflaton. 

Let us next consider the second term N(tf, t en d, 0"*)- When the slow-roll conditions are 
violated, the inflaton starts to oscillate around its minimum. After several oscillations, 
the energy density of the inflaton can be well approximated by that of a perfect fluid 
with a constant equation of state. If the potential around the minimum is written as 
V oc 2n , its energy density decreases as p^ oc a - 6ri /( n + 1 ). For the discussion of the 
curvature perturbations in the next section, we further decompose the second term into 
two parts. Taking t c as a time after several oscillations of the inflaton but well before the 
time of decay, N(tf, tend, 0"*) can be divided as 

N(t f , t cnd , a,) = N(t c , tend) + N(t f , t c , a*). (6) 

Note that the first term on the right-hand side does not depend on a* as long as t c is 
taken to be sufficiently before the time of the decay. As will become clear later, N(t c , tend) 
is irrelevant to the curvature perturbations. However, we have to take into account it 
when we calculate the epoch when the reference scale at present, which is taken to be 
k = 0.002 Mpc -1 for our analysis, exited the horizon during inflation. In general, we need 
a numerical calculation to evaluate N(t c , tend) and do not discuss this term further here. 
Regarding N(tf,t c ,a*), for our discussion in the following, we write it as 

N(t f , t c , cr,) = - log Bl + q [r^a*, t c )/H c ] , (7) 

4 p f 

where p c and pj are total energy densities of the universe at t c and t/ respectively. This 
equation should be understood as a definition of the function Q. Note that Q depends on 
the variable x = T(a^,t c )/H c alone, which can be confirmed by a dimensional analysis of 
the background evolution equations. Since the dependence of Q on x differs depending on 
the expansion law of the background space-time, i.e. the power 2n that determines the 
form of the inflaton potential, and also on the dominant decay channel of the inflaton, we 
need to follow the background evolution for each case. In principle, we have to evaluate it 
by numerical calculations. However, we can make approximate analytic estimates for Q(x) 
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Table 1: We list the analytic values of ao defined by Q(x) = ao logic for various infiaton 
potential V(<f)) oc <p 2n and dominant interactions. Notice that these values are in good 
agreement with numerically obtained ones, in particular for small values of x. If n = 1, 
the infiaton cannot decay only with the four-point interaction —h<j) 2 x 2 - 



when x< 1. We found that Q(x) can be well approximated with the form Q(x) = ao log x 
where ao is a numerical constant which depends on the interaction between the infiaton 
and matter and the potential for the infiaton. In Table (TJ the values of ao are listed for 
three different powers for the infiaton potential and three different interactions, which are 
obtained by analytic estimates where the sudden-decay approximation is adopted. The 
derivation of these analytic values is given in Appendix [B] We have also calculated them 
numerically and found that the differences between analytic estimate and numerical one 
are at most 10 % for x < 10~ 6 , which can be seen from Fig. [1] where the relative error for 
the case with the inflation potential being V oc cffi and three types of interaction listed in 
Tabled 

In the following analysis, we consider the cases with n — 1,2 and 3. Here we briefly 
discuss the tendency of Q for these cases in order. When n — 1, the infiaton potential can 
be written as V(<fr) = ^-0 2 . In this case, we have m e F = m and the decay rate becomes 
independent of time for the interactions such as —y(pijjijj and — M<px 2 - Thus we obtain the 
same value of Q for these interactions. Meanwhile, if the dominant decay occurs through 
the four-point interaction, the infiaton cannot decay completely into radiation simply 
because the decay rate which is proportional to the Hubble parameter squared decreases 
faster than the Hubble parameter. Thus in this case, the universe cannot be reheated. 
Hence we do not consider a four-point interaction case for the quadratic potential. When 
n = 2, P0 behaves the same as the energy density of radiation. Hence the universe 
expands in exactly the same way as the background space-time even if the decay rate of 
the infiaton slightly deviates from the background valua^H This means that Q = for 
any interactions. For the case with n = 3, since the effective mass rrfF changes in time, 
the decay rate evolves differently for different types of a dominant interaction, which gives 
a different form of Q. 

Also note that since the coefficient of logx in Q{x) is negative (positive) for n = 1 (n = 
3), the function Q(x) itself becomes positive (negative) for n — 1 (n = 3) when x < 1. 
From the definition of Q (see Eq. ([7])), Q represents the deviation of the e-folding number 



# 1 This is not true if preheating occurs. In Ref. [24], it was shown that the universe can evolve differently 
at different locations separated by the super-horizon distance during preheating because of the persisting 
isocurvature perturbations at the end of inflation. 



•a 



0.1 



0.01 



0.001 




1Q -20 10 -18 10 -16 10 -14 1Q -12 1Q -10 1() -. 



10"" 10"" 10"' 



Figure 1: The relative error for the function Q(x) between the one obtained by analytic 
and numerical methods for the cases with the interaction —ytpipip (red solid line), —M<px 2 
(green dashed line) and —h<p 2 x 2 (blue dotted line). Here we assumed V oc 6 for the 
inflaton potential. For the cases with the quadratic potential V oc 2 , the errors are 
smaller than those given in this figure. 



from the radiation dominated universe. Because energy density of dust decays more slowly 
than the radiation energy density, the more dust gives the more e-folding number until the 
total energy density decreases to a fixed value, which means the positive Q. Meanwhile, 
because the energy density of the inflaton oscillating in the sextic potential decays as a -9 / 2 , 
it decreases more rapidly than the energy density of radiation. Hence Q becomes negative 
for the sextic potential. 

In closing this section, we write down the expression for N(tf, £*, 0*, a*): 



N(t f ,U,(j)*,a*) 



Ml 



£ + ~ log Bl + AT(t C5 tend ) + Q \r+{ff<)/H ] , (8) 
v 4> 4 Pf 



which forms the basis in calculating various observables in the subsequent sections. 



4 Observables 



4.1 Curvature perturbation 

To calculate the curvature perturbation, we make use of the 5N formalism [20-23] . In this 
formalism, the primordial curvature perturbation ( on the uniform energy density hyper- 
surface at the time t — tf is given by differentiating the e-folding number N(tf, £*, 0*, a*) 



in Eq. (jSJ) with respect to 0* and a*, 



Ki v t> 2M p 2 i V v l J 6M p 2 i V v* v l v l 

+ QM* + ^QoMl + -QaaaSo-l (9) 

z o 

where <50* and 5cr* are the perturbations of 0* and cr* on the flat slicing. For the purpose of 
this paper, we include the terms up to cubic order in the perturbations of scalar fields. If 
we choose tf well after the reheating, then £(£/) gives primordial adiabatic perturbations. 
Using x = T </,(&*, t c )/H c , the derivatives of Q with respect to a can be expressed using 
r<£ as 

Q(7 =xQ'(x)^f = A(x)^f, (10) 

Qa* = xQ'(x)^f + x 2 g"(x)S = A(x)^f + B(x)S, (11) 

Q affa = xQ'(x)^ + Zx 2 Q\x) V -^p- + x 3 Q'"(x)^ 

= A(x)^ + 3B(x)^^ + C(z)G (12) 

Here we defined the function as A(x) = xQ(x)', B(x) = x 2 Q(x)" and C(x) = x 3 <5(x) w and 
denote the derivatives of T^ with respect to a* as T a , T aa , and r CTCrcr . When i< 1, which 
we consider in the followings, these functions become almost independent of x, namely Q 
can be well approximated by Q(x) = a logx (see Table [TJ). In this case, the functions 
defined above become constant and can be written as A = Oo, B = —a Q and C = 2a$. 

4.2 Power spectrum 

The power spectrum P^ of the curvature perturbations is defined by 

(Q i ( l : 2 ) c = (2rr) 3 P c (k 1 )5(k 1 + k 2 ), (13) 

where (■ • • ) c means that we take connected parts of (•■■)• 

By using Eq. (Q, we can express V^(k) = k 3 P^(k) /27r 2 to the leading order in 50*, 5a* 
as 
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p <<*>=*U^.)< 1+fl >- (14) 



where e is a slow-roll parameter defined by, 



^^fr US) 



Furthermore R = 2eA(x) 2 M^T 2 /T 2 is the square of the ratio of the curvature perturbation 
from modulated reheating to that from the inflaton, i.e., Cmod/Cilf- Thus the limit R — > 
(R —>■ oo) corresponds to the case where the curvature perturbation is sourced only by the 
inflaton (modulus) fluctuations. For discussion in the following, here we give the definition 
of another slow-roll parameter rj: 

V = M^. (16) 

4.3 Bispectrum 

The bispectrum B^ is defined by 

(CfcCfeCj^c = (2nfB c (h, k 2 , k 3 )S{k 1 + k 2 + k 3 ). (17) 

During inflation, both (p and a are slowly-rolling. It is known that the non-Gaussianity of 
£ coming from the intrinsic non-Gaussianities of 5(f)* and 5a* is far below the observational 
sensitivity. Hence we can treat 5(p* and 5a* as uncorrelated Gaussian random fields with 
the same amplitude. 

If we parameterize B^ by the dimensionless parameter /nl by 

B c (k u k 2 , h) = jj/ NL (P c (fc 1 )P c (fc 2 ) + P c (k 2 )P c (k 3 ) + P c (k 3 )P c (h)) , (18) 

then / NL can be written as 



6 r 2 ( b(x) rr 

5 /NL ~ A(x)(l + Rf \A(x) + Tl 



7/nl = -rTT7-^72 ( T77V + "FT ) + °^ V)- (19) 



The first term A(x)~ 2 B (x) R 2 / (1 + i?) which is independent of how the decay rate V 
depends on a represents the non-Gaussianity coming from the non-linearity between C, 
and 5Y. Because the function R 2 /{I + R) is suppressed by R 2 for R <^ 1 and approaches 
1 for R 3> 1, the magnitude of the first term is at most |v4~ 2 I?| = (9(1 ~ 10). For example, 
if the potential is given by a quadratic term, we have |t4~ 2 £>| ~ 6 (see Table [T]), which 
yields /nl = 5. The non-linearity between ( and 5T gives the positive (negative) /nl for 
quadratic (sextic) inflaton potential. 

On the other hand, the second term which represents the non-Gaussianity coming from 
the non-linearity between T and a depends on the detailed form of T(a). We see that very 
large non-Gaussianity |/nl| ~^> 1 can be generated only when irro-oyr 2 ! ^> 1 is satisfied. 

4.4 Trispectrum 

The trispectrum T^ is defined by 

(Cj^Cj^c = (2tt) 3 T c (A;i, h, k 3 , h)5(h + k 2 + k 3 + h). (20) 
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Figure 2: Consistency relations among three non-linearity parameters are shown for the 
case with modulated reheating corresponding to Eq. fl25l) (red solid line) and the curvaton 
model for the case where the curvaton decays before it dominates the universe correspond- 
ing to Eq. (77) in Ref [19] (black dashed line). In this figure, the consistency relation is 
presented as contours of #nl in the /nl _ Tnl plane. The line for models with "ungaussi- 
ton" [26] is also shown (blue dotted line). For this model, the relation between /nl and 
r NL is given, thus it is irrelevant to the value of #nl- Notice that the inequality Eq. (124)) 
should hold for the above mentioned scenarios. Thus we also show a region violating this 
inequality with shade. 



If we parameterize T^ by the two dimensionless parameters tnl and #nl as 

T((ki, fa, fa, fa) = t nl {P({kia)P((fa)P((fa) + 11 perms.) 
54 



+-£ NL (P c {fa)P^fa)P^fa) + 3 perms.) 



then these parameters are given by 



(21) 
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From Eqs. (fT9l) and (122|) . neglecting the terms of the slow- roll order, we find that 

^6 

r NL > |/^ L . (24) 

In Ref. [25], it is shown that this inequality holds in more general setting, i.e. for any 
cases where the leading non-Gaussianity comes from super-horizon evolution. 

As for (7nl 5 the non-linearity between ( and 5T always gives positive values for g^L 
irrespective of the inflaton potential and the dominant interaction. When the non-linearity 
between 5T and 5a cannot be neglected, #nl can be negative. 

If the T aaa term in #nl can be neglected and x <C 1, all the three non-linear parameters 
are functions of R and TT au /T 2 a . This means that we have a universal relation among the 
non-linearity parameters independent of R and TT au /T 2 a . The explicit form of such a 
relation can be written as 



864 f« L 12 /g 
625ag r^ L 5a r NL 



#NL - - cot .o 3 - — ^- l^5J 



In particular, if the potential around the minimum is quadratic, then we have a 
and the relation above reduces to 



31104 j% k + 72 fl 

625 t-^l 5 r N L 



^nl = -^^3^ + ^^. (26) 



By checking this consistency relation from observations, we can discriminate this scenario 
from other ones that also generate large non-Gaussianity. In Fig. El we show the contours 
of #nl as a function of /nl and Tnl given in Eq. (1261) . For comparison, we also show 
the relation among the non-linearity parameters which holds for other scenarios which 
can generate large non-Gaussianity such as models with mixed inflaton and curvaton for 
the case with a large positive /nl [19] and "ungaussiton" [26]. Interestingly, while the 
curvaton gives negative g^L, the modulated reheating gives positive one (if the potential is 
quadratic) if /nl is large positive. Hence just the determination of the sign of #nl enables 
us to discriminate these two competitive models in this case. We also plot the relation 
between /nl and tnl in the "ungaussiton" model where the relation between tnl and /nl 
is given irrespective of #nl- Thus only one line is drawn for this model. 

4.5 The scalar spectral index and the tensor-to-scalar ratio 

In this subsection, we will give the expressions for the scalar spectral index and tensor-to- 
scalar ratio. First of all, the scalar spectral index n s is given by 



lis 



d\ogV t 



k=aH 



dlogk 

4e - 2r? 

10 



During inflation, the tensor mode is also generated and its power spectrum can be 
written as 

H ^ 2 



To characterize the size of the tensor mode, the tensor-to-scalar ratio is usually used, which 
is given by 

r = ^ = ^ (29) 

5 Inflationary predictions and comparison with ob- 
servations 

In this section, we discuss the five inflationary parameters, the scalar spectral index, tensor- 
to-scalar ratio and three non-linearity parameters (n s ,r, / NL , tnl^nl) in models where 
fluctuations from modulated coupling and inflaton can both contribute to the primordial 
curvature perturbations. With regard to n s and r, current cosmological observations give 
severe constraints on these quantities, thus we also compare the predictions for n s and r 
with the WMAP 5-year data [7,27]. For the purpose of this paper, we consider the chaotic 
inflation models with some polynomials for definiteness. However, before we discuss its 
inflationary predictions, here we make some comments on the cases with other inflation 
models in a mixed scenario. If we take the new inflation models, its effects are expected 
to be very small since the new inflation models usually give a negligibly small value of the 
slow-roll parameter e. Since e appears in R = 2eA(x) 2 M^T 2 /T 2 which represents the size 
of the contribution from fluctuations from modulated reheating, R becomes very small 
when e is negligibly small. This fact has already been pointed out in Ref. [19] for a mixed 
model of the inflaton and curvaton. Hence, it is unlikely that fluctuations from modulated 
reheating dominate the total curvature perturbations in the new inflation model, in which 
we recover the usual standard formula for the inflationary parameters and the effects of 
modulated reheating can be negligible. As another possible inflation model, we can also 
assume the hybrid inflation model. However, as is also discussed in Ref. [19], the case of 
hybrid inflation leads to the similar result as that of chaotic inflation. Thus we omit them 
in this paper. Regarding the interactions between the inflaton and matter, we adopt the 
Yukawa interactions, just to be concrete. We can straightforwardly analyze the case with 
other interactions by simply replacing the numerical factor A(x), B(x) and C(x) appearing 
in the expression with the corresponding values (see Tabled]). 

Even after we fix the form of the inflaton potential and the interactions for reheating, 
we have still a degree of freedom of how the coupling constant y depends on the modulus 
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a. In the following, we assume that y a {cr) can be written aa^ 2 . 

y a (v) = y° a (i + a a ^ + p a ^ + •••), (30) 

where a a and (3 a are 0(1) coefficients. M is some energy scale and we assume that 
\a\ <C M. To be definite, we truncate the expansion of y a (a) at the second order in a. 
Then the decay rate can be also truncated at the quadratic order in a, 

r = r ,( 1 + a £. + p*L) l (3i) 

where a and (3 are also 0(1) coefficients. 
From Eq. (l3Tj) . we have 



VY aa 2(3 



n ' ( M 2(3a 



(32) 



If a is 0(1), then this equation further reduces to ~ 2 j3 /a 2 . From Eq. (|19j) . the non- 
linearity between T and a gives /nl — 2{3/(A(x)a 2 ) for _R ^> 1. Since |y4(x) _1 | is (9(1 ~ 
10), |/nl| = 0(10 ~ 100) can be achieved by setting 2f3/a 2 to be O(10), which is quite 
possible while satisfying a, f3 — Oil). Notice that, to obtain a large positive /nl, P should 
be negative for the quadratic potential since A(x) < for the potential. For example, if we 
take (a, (3) = (0.5, —1) and assume the quadratic potential and the Yukawa interaction for 
the inflaton in which A(x) ~ —1/6, then we have /nl = 45 for R 3> 1. Since the constraint 
on /nl from WMAP 5-year data is given as —9 < /nl < HI [7], this sort of possibilities 
may be interesting. Furthermore this demonstrates that the non-linearity between T and 
a can provide non-Gaussianity of £ much larger than that from the non-linearity between 
£ and ST. 

Meanwhile, some symmetries may forbid the appearance of the linear terms in a in 
ya{c)- In this case, Eq. (132]) becomes ~ M 2 /(2/3<t 2 ), where /nl can be very large because 
of M/a ^> 1 when /3 < 0. For an illustrational purpose, we will consider two cases: 
(a, /3) = (0.3, -1.0) (case A) and {a, 13) = (0.0, -1.0) (case B) in the following. 

To compare the prediction for the primordial curvature fluctuations and non-Gaussianity 
with observations, we need to specify when the present cosmological scale exited the hori- 
zon during inflation. Since fe* = a{t*)H* holds when the scale with the wave number fc* 
crossed the horizon, the reference scale fc re f where we probe the primordial fluctuations at 
the present time is related to that at the horizon crossing during inflation as 

&ref a{U)H* 



a H a H 



(33) 



# 2 Here we assume that there are no renormalizable interactions between a and ip a , Xa because such in- 
teractions will give the modulus a thermal mass much larger than the Hubble parameter, which suppresses 
fluctuations of modulus and spoils the modulated reheating scenario. 
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where a^ and Hq are the scale factor and the Hubble parameter at present. The ratio of 
a* to a on the right-hand side can be divided into several parts as, 

fcref a{t*) a(W) a(tf) H* 

O-oHq O (tend) «(*/) a H 

The definition of t ea & and tf is given in section [3j By taking logarithm of both sides 
and using Eqs. © and (J7|), the number of e-folding between the time t* and t e nd, i-e., 
N(t en d, £*, 0*) can be written as 

iV(W, t*, 0*) = - log -%- - N(t c , t ead ) - - log ^ - Q(a;) + log ^ + log ^, (35) 

with which we can determine the field value 0* at the horizon crossing. 

For the reference scale, we take k rc { = 0.002 Mpc -1 in the analysis. As for the fifth 
term, we assume that no more entropy is produced after the inflaton decays. Thus this 
term can be rewritten by using the conservation of the entropy density per comoving 
volume. Since the entropy density is given by s = (2ir/4:5)g ifS T 3 with g* s being the total 
number of effective massless degrees of freedom, we have the following relation, 

log— = -log— = -log -£. (36) 

a 3 s f 3 g* s fTf 

For gt, s f at the time of a/, we take #* s / = 100. With regard to other quantities, we assume 
H(t c ) = 10~ 2 H cn d and x = 10 -8 in the following analyses. 



,2 



,2 



5.1 Chaotic inflation : V{4>) = ^ ( 

Now let us first consider chaotic inflation [28] with the quadratic potential^ 3 !. Since we 
have fixed the values of the quantities which determine the background evolution after 
the end of inflation as mentioned above, remaining variables that we need to specify 
are some parameters in the inflaton potential and the values of a and M relevant to 
fluctuations from modulated reheating, which appear in the decay rate of the inflaton. 
Since the inflation is assumed to be driven solely by the inflaton, the Hubble parameter 
during inflation is controlled by the value of the parameters in the potential. Thus the 
primordial curvature fluctuations also depend on these parameters, which means that the 
parameters in the inflaton potential can be fixed by the WMAP normalization, i.e. by 
requiring that the total curvature fluctuations are O(10~ 5 ). In fact, the normalization 
slightly depends on the spectral index and tensor-to-scalar ratio. Hence we used the one 
given in Ref. [36j^_jderived using the WMAP 5-year data in the following analyses to fix a 



# 3 This type of a simple polynomial potential can be realized in supergravity [29-35] . 
# A The amplitude at k = 0.002 Mpc" 1 is given with 8 2 H = (4/25)P c as 

105 ,wmap 5 = L91Q x cxp [(-0.724 + 0.533 r)(l - n,)] 

H VI + 0.278 r K ' 

For the details, see Ref. [36]. 
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Figure 3: Contours of R in the a-M plane for the chaotic inflation model with quadratic 
potential. Left and right panels are for case A [a = 0.3, /3 = —1) and case B (a = 0,/3 = 
—1), respectively. In the figure, M and a are shown in units of M v \. 



parameter in the potential for the inflaton. For the quadratic case of the chaotic inflation, 
we fix the value of m by the WMAP normalization. With regard to the spectral index n s 
and tensor-to-scalar ratio r, the observation of WMAP 5-year also give severe constraints 
on these quantities (see Fig. 3 in Ref. [7]). When fluctuations from the inflaton alone are 
responsible for the curvature fluctuations, in our setting (fixing the parameters) where the 
number of e-folding during inflation is A^ n f ~ 60, these quantities are given as n s = 0.967 
and r = 0.133 which are allowed by WMAP 5-year data. However, if fluctuations from 
the modulus fluctuations are included, these predictions can be modified, which we will 
discuss in the following. To see in what cases such modifications are significant, we show 
the contours of R in the a-M plane in Fig. [3]for the case A (left panel) and B (right panel). 
By looking at the figure, we can expect the parameter region where the predictions for n s 
and r are modified much. 

In the top panels of Fig. IH we show contours of the scalar spectral index n s in the a-M 
plane for the case A (left panel) and B (right panel). Let us first look at the left panel 
(case A). In this case, Y a /Y is almost independent of a . Hence R is also independent of 
a and depends only on M. As can be seen from Eq. (127)) . the effect of fluctuations from 
modulated reheating on n s appears only through the parameter R. This is the reason why 
the contours are parallel to the a-axis. Furthermore, since R ex M~ 2 , smaller value of 
M indicates larger R, i.e. larger contribution from the modulus fluctuations to the total 
curvature perturbations. When the modulus contributions become dominant, the value of 
n s shifts from 0.967 to 0.983. For the case B (see the top right panel of Fig. HJ), T a /T is 
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Figure 4: Contours of n s (top panels) and r (bottom panels) in the a-M plane for the 
chaotic inflation model with quadratic potential. Left (right) panels are for case A (B). 



proportional to a/M 2 . Hence the contours are parallel to the line M/a 1 ^ 2 = const., which 
can be seen from the figure. 

Now we investigate how the tensor-to-scalar ratio is modified in the mixed model. In 
the bottom panels of Fig. HI we show contours of the tensor-to-scalar ratio r in the a-M 
plane. Notice that as in the case of n s , the effect of the modulus on r appears only through 
R. Hence the slope of the contours becomes the same as those for n s . We see that when 
the modulus contributions are dominant, r becomes negligibly small. 

Next we move on to the issues of non-Gaussianity. For this purpose, we show contours 
of non-linearity parameters /nl 5 tnl and #nl i n the top, middle and bottom panels of 
Fig. El respectively. Remember that, when the fluctuations from the inflaton alone are 
assumed, these non-linearity parameters are of the orders of the slow-roll parameters, which 
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Figure 5: Contours of /nl in the a-M plane for the chaotic inflation model with quadratic 
potential. Left (right) panels are for case A (B). 
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are usually much less than unity. However, since fluctuations from modulated reheating 
can give large non-Gaussianity, the non-linearity parameters can be large even when the 
contribution from the modulus fluctuations is subdominant in the curvature fluctuations. 
For the case A, we have TT aa /T 2 ~ const, (see Eq. (1321 ). Hence as in the cases for n s 
and r, the value of R alone determines the effects of the modulus contributions on the 
non-linearity parameters. Thus the smaller M gives the larger non-Gaussianity which 
can be seen from the figures. On the other hand, for the case B, the above mentioned 
combination is TT aiJ /r 2 ~ M 2 /a 2 . Hence a and R can both affect the non-linearity 
parameters in this case. If the former effect is dominant, the slope of the contour becomes 
1/2. On the other hand, if the latter one is dominant, its slope becomes 1 and we see 
that for larger R and smaller M 2 /a 2 , we have larger non-Gaussianity. In fact, in Fig. [5j 
contours corresponding to the latter case are not shown for this model, since the size of 
the non-linearity parameters become too large in such a region for this model, thus we do 
not depict them here. However this kind of behavior can be seen in Figs. fTUl and fT2l which 
are for other inflation models to be discussed in the following. 

For an illustrational purpose, we fix the values of a and j3. However, the predictions 
for the inflationary parameters also depends on these values. Thus here we stop to discuss 
its dependence. In Fig. [61 contours of the ratio R are shown in the a-M plane. Since the 
relative size between a and (3 is important to see the effects of the modulated reheating, 
here we fix the value of (3 and vary a. Since the ratio R is controlled by the combination 
T^/r 2 , as a increases, R also becomes large for M being fixed. As discussed in the previous 
section, n s and r are also determined by the above combination, thus the tendencies are 
the same as that of the ratio R, which can be seen in Fig. [7] where contours of n s and r 
are depicted. 

However, if we look at plots of non-linearity parameters, which are shown in Fig. [8], 
the trends are different. As mentioned before, the non-linearity parameters are governed 
by the combination TT aa /T 2 a . ~ 2(3 /(a + 2/3a/M) 2 , when a is small, the dependence of 
non-linearity parameters on a is also small. On the other hand, when a is large, the size 
of the above combination is determined by a with (3 being fixed. Thus the predictions 
are irrelevant to M in this case, which can be read off from the figure. In the following, 
we discuss other types of the inflaton potential and again present our results in the a—M 
plane fixing the values of a . However, the trend discussed here also applies to those cases. 

5.2 Chaotic inflation : V{(f>) = f 4 + ^ (\? 

Next we consider chaotic inflation with the quartic potential. As a simple possibility, one 
usually assumes that the quartic term alone exists in the potential for the inflaton. Since 
the energy density of inflaton oscillations in such a potential decreases in the same way as 
that of radiation, no fluctuation is generated via modulated reheating scenario as discussed 
previously. However, there are some interesting models [37,38] where the potential during 
inflation is given by the quartic form, but the quadratic term becomes effective when 
the inflaton starts to oscillate. In fact, this feature comes from the curved trajectory 
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Figure 6: Contours of R in the a-M plane for the chaotic inflation model with a quadratic 
potential. 
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chaotic inflation model with a quadratic potential. 



in the multiple fields configuration. But it can be well characterized by considering a 
single field with the potential given by the sum of the quadratic term and the quartic one: 
V(4>) = ^-0 2 + 1 4> 4 - Above situation can be realized by taking m = \/\M p \ with A being 
fixed by WMAP normalization. Hence we adopt this relation in the following analysis 
for definiteness. In this case, quartic term dominates during inflation (0 > M p i), but the 
quadratic one dominates after inflation. 

This model is also interesting in another aspect. In the absence of the quadratic term, 
we have Q = and the curvature perturbations are not generated by the modulation of the 
decay rate. Then the curvature perturbations originate solely from the inflaton fluctuations 
even when the decay rate fluctuates, which is completely ruled out from WMAP 5-year 
data because of the too large tensor-to-scalar ratio. However, if we include the quadratic 
one, the curvature perturbations can be additionally generated when inflaton decays. If 
these perturbations dominate the total curvature perturbations, the quartic inflation model 
may still satisfy the observational constraints, with large amount of the non-Gaussianity 
which could be tested by the future observations. 

Fig. [9] show contour plots of n s and r in this case, respectively. Shaded regions are 
excluded by the WMAP 5-year constraints on (n s ,r). We see that when the modulus 
contributions become dominant, n s gets closer to unity and r becomes very small. As a 
result, (n s , r) enters the allowed region of WMAP 5-year data. The slopes of the contours 
are essentially the same as in the previous case where V(4>) = ^r4> 2 - 

In Fig. [TU1 we show the contour plots of the three non-linearity parameters. As for the 
case A, there is a large parameter space where (n s ,r) is in the allowed region of WMAP 
5-year data while generating large amount of non-Gaussianity /nl = O(10 ~ 100) which 
could be tested by the upcoming observations. 
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Figure 9: Contours of n s (top panels) and r (bottom panels) in the a— M plane for the 
chaotic inflation model with quadratic and quartic potential. Left (right) panels are for 
case A (B). Shaded region is excluded by WMAP5. 



5.3 Chaotic inflation : V(<f>) = f(^) 

Let us next consider chaotic inflation model whose potential is given by 

v(0) v " 



6 V*W ' {X 

This model is also completely ruled out by WMAP 5- year data if the curvature perturba- 
tions solely originate from the inflaton fluctuations. 

Fig. [TT] show the contour plots of n s and r in this case, respectively. We see that as in 
the case of the quadratic and quartic potential, (n s , r) enters the allowed region of WMAP 
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the a-M plane for the chaotic inflation model with quadratic and quartic potential. Left 
(right) panels are for case A (B). Shaded region is excluded by WMAP5. 
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Figure 11: Contours of n s (top panels) and r (bottom panels) in the a-M plane for the 
chaotic inflation model with sextic potential. Left (right) panels are for case A (B). Shaded 
region is excluded by WMAP5. 



5- year data due to the slight shift of n s to unity and the significant suppression of r when 
the modulus contributions become dominant. 

In Fig. [121 we show the contour plots of the three non-linearity parameters. For 
the sextic potential, /nl becomes negatively large when the modulus contributions are 
dominant, as discussed in the previous section. If we take f3 < instead of (3 > 0, then 
we have positively large /nl just like the case in the quadratic and quartic potential. 
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6 Summary 

We have investigated a mixed scenario where both fluctuations from the inflaton and 
modulated reheating are responsible for cosmic density fluctuations today. First we sum- 
marized the decay rate of the inflaton for various interactions and inflaton potentials, 
whose details are presented in Appendix |A] Then we gave expressions for the spectral 
index, tensor-to-scalar ratio and non-linearity parameters such as /nl, tnl and #nl- We 
found that while the non-linearity between ( and 5T generates /nl = 0(1 ~ 10), the non- 
linearity between T and the modulus can generate large /nl > 0(100). We also derived 
the consistency relation Eq. (125]) among the non-linearity parameters which is applicable 
for the cases where the decay rate for the inflaton is given as Eq. (I3TT) . If the potential 
around the minimum is quadratic, then $nl takes the same sign as /nl- This is in sharp 
contrast to the situation where the mixed model of the inflaton and the curvaton yields a 
negative #nl when /nl is large positive. 

We have also studied the inflationary parameters including non-linearity parameters 
assuming some inflation models. Since the new inflation model usually gives a extremely 
small value of e, we have R <C 1. Thus, the addition of fluctuations from modulated 
reheating does not significantly change the prediction of the property of fluctuations in the 
case of new inflation. Therefore, in this paper, we have considered chaotic inflation models 
with several polynomials and some types of decay rate for an illustrational purpose. In 
the case with the chaotic inflation model, the contribution from fluctuations of modulated 
reheating makes non-Gaussianity larger, the spectral index closer to scale-invariant and 
tensor-to-scalar ratio more suppressed. This helps to liberate models of inflation such as 
the sextic potential one because this model has been excluded by the data because of its 
too red-tilted spectrum and too large tensor-to-scalar ratio. The chaotic inflation with 
quartic potential has also in fact been excluded by the data. Since fluctuations are not 
generated by modulated reheating in a simple quartic potential case, this model cannot 
be made viable by just adding a contribution from the modulated reheating. However, 
even if the inflaton is driven by a quartic term, another term such as a quadratic one 
can become effective during reheating stage after inflation. This kind of situation can be 
realized by assuming V(<p) = (A/4)0 4 + (l/2)m 2 </> 2 with m = vAM p i, which was discussed 
in this paper. In such a case, fluctuations from modulated reheating can affect the total 
fluctuations and liberate the model. In addition, non-Gaussianity can also be large in 
this case too. We have also investigated the case with the quadratic inflation model and 
showed that non-Gaussianity can be large in this model by just adding some contribution 
from modulated reheating without conflicting the observations of the spectral index and 
tensor modes. 

In most works so far, observational consequences for various generation mechanisms of 
primordial fluctuations have been discussed in each separately. However, different mecha- 
nisms can be in general responsible for today's cosmic density fluctuations simultaneously 
so it is of great importance to investigate a mixed model such as the one we discussed 
in this paper. Precise cosmological observations expected in the near future to hunt for 
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primordial non-Gaussianity may reveal such interesting scenarios in the early universe. 
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Appendix 

A Calculation of the decay rate 

In this appendix, we calculate the decay rate of the inflaton to lighter particles for three 
different types of interactions [39-41]. We also consider three different types of potentials 
around the minimum: quadratic, quartic and sextic potentials. In the following, we ap- 
proximate the background space-time as Minkowski space because the energies of created 
particles from the inflaton are much higher than the Hubble parameter. In this appendix, 
we follow the notations of Ref . [42] . 

A.l Yukawa interaction: C\ nt = —y(fnjjijj 

We regard <p as the classical field which is spatially homogeneous but oscillates in time t. 
Then the interaction part of the Hamiltonian is given by 

V = y</>(t) J d 3 x^. (39) 

ip in the interaction picture is equal to that in the Heisenberg picture when the interaction 
is absent. Hence we write if) as 

^ = E / T^b ( tt ® °") e * a v + v & ^ e ~ ipxb t) > ( 4 °) 

a J ( 27r ) V ' 

where a denotes the helicity and bf >cr and b c =* a satisfy the following relations, 

Ofc*. t>U = {6£„, b c l,} = S(p- g)<W, {6 P> , b^} = {*&,, blA = 0. (41) 
Also, u(p,a) and v(p,a) satisfy the following equations, 

{ip^ln + m^)u(p, a) = 0, (-ip^,, + m^)v(p, a) = 0. (42) 

Let us write the period of the oscillations of <f>(t) as T and Fourier expand cj)(t) as 

oo 

<P(t)= E ^e~ lwn \ (43) 

n=— oo 

where uj = 2%/T. 

Using these, the transition amplitude from the initial state \i) = |0) at t — — oo to the 
final two-particle state |/) = btb c J,\0) at t = oo is given by 



/oo °° 
dt(f\V(t)\i) = -2my5(p + q)u(p,a)v(-p,a') J^ <f> n 6(2E p - 
-oo __ , 



nu), (44) 



26 



where E p = Jp 1 + m, is the energy of the ip particle. Then T, the transition rate per 
unit time and unit volume, is given by 

T = V 2 y2\ ( t ) n\ 2 / — — 1 5{2E p -nu) y2\u(p,a)v(-p,a')\ 2 

n=l J ( 27r ) 

2 °° 

y 2 v^ 21 j. 12 

= *>>. (45) 

Here (■ ■ • ) denotes the average of • • • over one period of the oscillations. From the first 
line to the second one, we have used the relation, 

£|u(p>M-p>')| 2 = ^, (46) 

and assumed -0 is massless, i.e., m^ = 0. 

T is the production rate of two ^-particle from the vacuum, not the decay rate of 
0- field. The decay rate of the 4> energy, which we denote as 1^,, can be obtained from the 
energy conservation, 

p+T+At = ETAt. (47) 

The left-hand side denotes the energy loss of the 0-field during the infinitesimal time At 
while the right-hand side the energy gain of two ip-psaticles. E is the expectation value of 
the energy of the final two-particle state. From this equation, T^ can be written as 

r> = f E ^l. (48) 

87T p^ 

From the definition of E, it can be written as 



En,aa'I d3 P d3 lKP+q)^E f -na;)\M n \ 



E = ^\ n lr,~J e ,l TTTT^i (49) 



where Ef is the energy of the final state and Ai n is defined by 

M n = -27ryiu(p, cr)v(-p, a')(p n . (50) 

Substituting this into the equation above yields 

y^°° 1 \<t>r>\ (nu) , . 

E = ^ n = l |yn ' V L. (51) 

„=il0n| (nw) 
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Let us define the numerical factor a by 



T°° \6 IV 

= l^n=l Ignl » / 52 N 

n=l l^nl n 2 



Then T^, can be written as 



T^ = — u;a . (53) 

Hence once we specify the motion of (f)(t), we can immediately calculate T^ using Eq. ( [531) . 

A. 1.1 Quadratic potential:^) = ^-<p 2 

In this case, </>(£) can be written as 

(f>(t) = (p cos(mt), (54) 

where the frequency, m, is just equal to the mass of the inflaton. 
Then, we have 

a = 1, — = 1. (55) 



Hence T^ is given by 



A. 1.2 Quartic potential:l/(0) = 

In this case, <f>(t) can be written as 



on 



¥' 



where the frequency u; is given by 

a; = l/Z£® m f. (58) 

2V6r(|) v ; 

m^ ff is the effective mass of defined by 

mf = y/3Xc/) . (59) 
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Using these, we find 



3 
Then the decay rate of the inflaton can be written as 

V 2 _,eff 



a « 1.036, ^-^ = -. (60) 



r = A 2 f-m e /, (61) 



where A 2 ~ 0.676 is a numerical constant. 

A. 1.3 Sextic potential:^) = ^4> 6 

In this case, the frequency of the inflaton oscillations is given by 



-=1S<' (62 > 



mf is the effective mass of defined by 

m * = ~P r ^ ' (63) 

To determine a, we numerically solved the equation of motion for <fi over one period of 
the oscillations. The result is 

(</> 2 > 3 , N 

a « 1.0897, ^— ^ = -. (64) 

P</> 2 

Then the decay rate of the inflaton can be written as 

I> = A 3 fmf, (65) 

where A3 m 0.546 is a numerical constant. 

A. 2 Interaction with scalar field \ : £int = ~M(f>xx 
In this case, the interaction part of the Hamiltonian is given by 

f = M<p{t) / d 3 x XX- (66) 



Let us expand x a s 

X{x) = I - ^jf 1 = [fap+ e-^4) . (67) 
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f 2 P ; U^ap+e-^al 

(27r) 3/2 v^ V 



Then the transition amplitude from the initial state \i) = |0) at t — — oo to the final 

.tfltl 
p q' 



two-particle state |/) = atat|0) is given by 



/co (h 

dt (f\V(t)\i) = -2mM 5(p + q) ^ ^-S(2E p - nu). (68) 



The transition rate per unit time and unit volume becomes 

M 2 °° 

re=l 



57T 

The corresponding decay rate of the inflaton is given by 



^T<^>- ( 69 ) 



r, = f *£> , (70 ) 

8tt p^ 



where £?, the mean energy of the two-particle state, is given by 



Introducing the dimensionless number j3 by 

Eoo 
n=i n 

££=1 |0n 



2 



£ = ^n=i I^IL. (71) 

Eoo I / \2 V / 

n=l IP«I 



/^%^Trl. (72) 



r^ becomes 



r, = £*,<#. (73) 

8tt p 



A. 2.1 Quadratic potentials ((/)) = \<\? 
Using Eq. (1531) . we have 



p<f> m 
Hence T^ is given by 



— = i- (74) 



M 2 

r* = E:- (75) 



57rm 
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A. 2. 2 Quartic potential:V(0) = f 4 

Using Eq. (137)1 . we find 



P<t> (mf) 
Then the decay rate of the inflaton can be written as 



2 ) 5 48 
(3 ~ 1.004, ^-l « --^. (76) 



57rm 
where I?2 ~ 2.693 is a numerical constant. 



M 2 
r </> = ^2^-— if, (77) 



A. 2. 3 Sextic potential^ (0) = 

In this case, we numerically found 



6M*' 



P4> (mf) 
Then the decay rate of the inflaton can be written as 



2\ 12 Q3 
0« 1.010, —^—^. (78) 



M 2 
r = B,—-^ , (79) 



$7rm3 ff 



where .Bq « 4.362 is a numerical constant. 



&*,2 



A. 3 Interaction with scalar field %: £j nt = —h(f>x^ 
In this case, the interaction part of the Hamiltonian is given by 

V = h<p 2 (t) Jd 3 xxx- (80) 

Let us define ((t) as 



! (*) - <^ 2 > = E ^ e_ 



incut 



Note that u; is the frequency of C(£). Then the transition rate per unit time and unit 
volume is 

n=l 

= ^ ((0 4 ) " (0 2 ) 2 ) • (82) 
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Introducing the dimensionless quantity 7 as 

•,00 i/. i2 



Z^n=l Isn 

the inflaton decay rate is given by 



7 = ^ ' ,i , (83) 



h 2 (0 4 )-(0 2 ) 2 , , 

T, = :~7 ^ V ■ (84) 



A.3.1 Quadratic potential:^) = \ 

Using Eq. (151]) . we have 



7 = 1, — WV/_ = J?o (g5) 



p^ 4m 

Hence T^ is given by 



32irm 167rm 3 

A.3.2 Quartic potential:V(» = |0 4 

Using Eq. (137)1 . we find 



h 2 (t? h 2 

^ = ^ = J—,Pf (86) 



(04\_/^2\2 05Q 

7 « 1.007, — ^- « — — . (87) 



A 
Then the decay rate of the inflaton can be written as 

h 2 
87r(mS , ) 
where C2 ~ 8.86 is a numerical constant. 



A.3.3 Sextic potential:V(0) = 

In this case, we numerically found 



YoAfi 



6M 



pi 



7W1 .019, ^W!»S (89) 



P4> Vo<Po 

Then the decay rate of the inflaton can be written as 



*V = ^3—^—3^0, (90) 



fF\3 "t>i 



where C3 ~ 37.26 is a numerical constant. 
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B Calculation of Q(x 

In this section, we will provide explicit form of Q(x) for three types of inflaton potential, 
quadratic V oc 2 , quartic V oc 4 and sextic V oc 6 , and for three types of interactions 
between the inflaton and matter particles, Yukawa interactions, three-point interactions 
£i nt = —M(fix 2 and four-point interactions £; nt = —hcj) 2 x 2 - 

B.l Quadratic potential 

If the inflaton potential is quadratic, i.e. V((f>) = ^-<fi 2 , then we have m cS = m. Hence 
the decay rate becomes independent of time for Yukawa interactions and three-point in- 
teractions £i n t = — M0x 2 . In this case, it was shown in Ref. [25] that Q(x) = — |logx. 
Meanwhile, if the dominant decay occurs through the four-point interactions, then the 
inflaton decay rate decreases as T^ oc p$. From the Friedmann equation H 2 oc p^, the 
ratio T^/H decreases in proportional to a _3//2 . Since T^ is smaller than H at initial time, 
that is, at the end of inflation, T^ never becomes larger than H. Hence the universe is 
never reheated only by the four-point interactions, which cannot realize the hot big bang 
cosmology. Thus we do not consider this case. 

B.2 Quartic potential 

If the inflaton potential is quartic, p^ decreases in proportional to a -4 . Hence the universe 
expands in the same way as the radiation dominated universe. From the definition of Q(x) 
(see Eq. (j7j)), we find 

Q(x) = 0, (91) 

for all three types of interactions. Note that unlike in the case of the quadratic potential, 
the universe can be reheated only by the four-point interactions in this case. 

B.3 Sextic potential 

If the inflaton potential is sextic, p^ decays in proportional to e~ 9N ^ 2 . Then the decay rate 
of the inflaton evolves as oc e~ 37V//2 , e 37V//2 and constant for Yukawa interactions, three- 
point interactions C mt = —Mifix 2 an d four-point interactions C mt = —h(j) 2 x 2 , respectively. 
Since the Hubble parameter decays faster than e~ 2N , the ratio T^/ H grows in time for any 
interactions. Hence the universe eventually becomes the radiation dominated universe. 

Since the most of radiation is produced when H = T^, which we checked by a numerical 
calculation, in what follows, we will use the so-called sudden decay approximation, where 
the inflaton decays instantaneously when the decay rate becomes equal to the Hubble 
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parameter. Under this approximation, the background equations can be written as 

%fi + 4pr = ^P+W-N d ), (92) 

r * P4> 5(N-N d ), (93) 



dN 2 rv H 1 

H 2 = ^r(p<p + Pr). (94) 

Integrating Eq. ( |92l) from N = which corresponds to t = t c to N = Nf > Nd, we find 

Wt, tc) = ^log P - + ^log ^ + 7V d . (95) 

Comparing this equation with Eq. ([7]), Q(x) can be written as 

Q = llo g T -^ + N d . (96) 

B.3.1 Yukawa interactions 

In this case, the decay rate can be written as T^ = r^(£ c )e _3Ar//2 . Hence from the equation 
H = r^,, we find that Nd is given by 

Nd= A iogI W- (9r) 

Substituting this into Eq. ( 1961) gives 

Q = ^log^. (98) 

o ti c 

B.3.2 Three-point interactions £ int = — M0x 2 

In this case, the decay rate is T^ = T ( j } (t c )e 3N ^ 2 . Then the corresponding Q is given by 

g = ll 0g £iM. (99) 

^ 30 & H c K J 

B.3.3 Four-point interactions C mt = —h(j) 2 x 2 

In this case, the decay rate is T^ = T^(t c ). Then the corresponding Q is given by 

Q = T% logT - j ir- (100) 
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